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1 Introduction 



One motivation for the recent resurgence of interest in massive gravity has been the need 
to explain the observed cosmic acceleration. Massive gravity, as an alternative to general 
relativity at cosmological scales, could provide the means to address the old cosmological 
constant problem (why the cosmological constant is zero) as well as an alternative explanation 
of the new dark energy problem (see [1] for a review). Here we investigate the possibility of 
a 'Partially Massless' theory gravity which can offer a new approach to the old cosmological 
constant problem. 

A generic theory of massive gravity requires a fiducial or reference metric upon which 
the theory of gravitons are defined. Maximally symmetric reference metrics are of special 
importance, as the notion of a graviton can then be defined via the irreducible representations 
of the associated isometry group. In Minkowski space-time, the representation theory of the 
Poincare group tells us that a massive spin-two field has five degrees of freedom. Generically, 
massive spin-two representations of the de Sitter group also have five degrees of freedom. 
However, for de Sitter there is also an exceptional case. For a specific value of the graviton 
mass relative to the de Sitter curvature, representation theory tells us that a massive graviton 
has only four propagating degrees of freedom. The field theory that describes this particle 
is known to exist at the linear level and is referred to as 'Partially Massless' (PM) gravity 
[2-9]. 

In the linear PM theory there is a scalar gauge symmetry which renders the helicity-0 
mode pure gauge and hence unphysical. If PM gravity exists as a non-linear theory (whether 
this is possible is the topic of this article), there will exist a single gauge symmetry which 
eliminates the helicity-0 mode at all orders. As with diffeomorphism invariance in general 
relativity, the PM gauge symmetry should be sufficient to fix the low energy form of the 
theory. This is in contrast to generic massive gravity for which no known symmetry enforces 
the form of the interactions. Thus PM gravity would have many of the aesthetic virtues of 
Einstein's theory despite being fundamentally an infrared modified theory of gravity. 

However, PM gravity differs from both general relativity and generic massive gravity 
in one crucial aspect. In PM gravity the value of the cosmological constant is not free, but 
is tied to the mass of the graviton via the gauge symmetry. To be precise, in D space-time 
dimensions, the cosmological constant A is given in terms of the graviton mass m by A = 
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{D — l)m? /2. The existence of a gauge symmetry ensures (in the absence of anomahes) that 
quantum corrections preserve this relation. This means that the old cosmological constant 
problem, which arises due to the large sensitivity of the value of A to quantum corrections, is 
replaced with the significantly more manageable problem of looking at quantum corrections 
to the mass of the graviton. The latter problem is under much better control due to the fact 
that in the limit m — )■ PM gravity reduces to general relativity which has an additional 
three gauge symmetries. The existence of these gauge symmetries in the massless limit is 
sufficient, via the t'Hooft naturalness argument, to protect the value of the graviton mass 
from large quantum corrections [10, 11]. The implication is that, if PM gravity exists as a 
non-linear theory, then it can provide a technically natural solution to the old cosmological 
constant problem which is based in symmetry and does not rely on degravitation or screening 
[12, 13]. 

PM gravity, if it exists, has a number of additional virtues over generic massive gravity 
that are a consequence of the absence of the helicity zero mode. There is no vDVZ discon- 
tinuity [14, 15] in the limit m — > provided that we maintain the PM relations between the 
mass and cosmological constant. In this theory there is no need for the Vainstein mecha- 
nism [16] - the extra scalar can't cause any fifth force problem because it simply doesn't 
exist - and hence no associated strong coupling. The theory would have a higher cutoff 
A2 ~ (Mpim)^/^ than that of generic massive gravity A3 ~ (Mpim^)^^^, as the worst non- 
renormalizable operators are associated with the now absent scalar field couplings. Finally, 
the absence of a hehcity zero mode could potentially remove any issues to do with superlu- 
minalities [17-19]. This is because superluminalities typically arise when the Galileon-like 
derivative couplings of the scalar are expanded around non-trivial backgrounds. In the PM 
theory these couplings would be absent. Combined with the potential of PM gravity to solve 
the old cosmological constant problem, these reasons compel us to look for the existence of 
the non-linear PM theory. (Ref. [20], which appeared as we were writing, asks this same 
existence question, and reaches the same conclusion we do, using different methods.) 

In the rest of this section we review the linear partially massless theory. We also review 
the ghost-free dRGT theories of massive gravity among which we search for a non-linear 
completion of PM gravity. We present a compact formulation of our candidate non-linear 
PM action, as was originally determined by the analysis of [21]. 
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In Section 2 we give independent evidence in support of our candidate action being the 
unique potentially partially massless theory. In particular, we perform an analysis for the 
vector modes in the decoupling limit, an analysis away from the decoupling limit using an 
FRW ansatz for the dynamical metric and its perturbations, and a brute force perturbative 
analysis of the non-linear PM symmetry up to cubic order in the fields. We find that the 
helicity-zero mode is indeed absent for the unique choice of the dRGT coefficients in the PM 
candidate theory. In particular, we give a detailed derivation of the following (known) result: 
while at quadratic order in fields the PM theory exists in arbitrary space-time dimensions, 
at cubic order it exists in only for D = 4 [22, 23]. We can also determine the non-linear 
gauge symmetry of the partially massless theory, in the FRW ansatz and to cubic order in 
interactions for a general metric. 

In Section 3 we show that, despite the promising evidence of the previous section, 
the candidate partially massless theory does not have the required gauge invariance to all 
orders. We see the first indication of this by looking at anisotropic cosmologies. We then 
push the perturbative analysis of the non-linear PM symmetry up to quartic order in the 
fields, making no assumption as to the form of the potential for the graviton. We show that 
at quartic order no PM symmetry exists without the introduction of either non-canonical 
derivative terms or additional fields. These conclusions agree with those of [20, 22]. The 
Lagrangian that most nearly realizes the PM symmetry at the non-linear level coincides with 
an expansion of the dRGT ghost-free Lagrangian with parameters given by the decoupling 
limit of [21], i.e., our candidate theory. Nevertheless, away from the decoupling limit an 
obstruction unavoidably remains. We conclude with a discussion of possible ways out and 
future directions for partially massless gravity. 

Conventions: We use the mostly plus metric signature convention, r]^^ = (— , +, +,+,.. .) 
and we work in arbitrary D > 3 space-time dimensions, unless otherwise stated. 
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1.1 The linear PM theory 



Consider the Fierz-Pauli theory of a massive graviton /i^;^ propagating on a maximally sym- 
metric background g^^, 

1. 



S 



R 



XI) 



Here th? is the graviton mass. The metric, covariant derivatives and constant curvature 
R are those of the background. For most choices of m? this action propagates the usual 
number of degrees of freedom of a massive graviton (e.g., five for D = A). For iv? = there 
is linearized diffeomorphism symmetry 6h^^ = + ^u^fi, and the action propagates the 
degrees of freedom of a massless graviton (e.g., two for D = 4). 

A clean way to see the degrees of freedom is to introduce, following [10], the Stiikelberg 
fields Afj^ and through the replacement 



"Ill/ 



/J,-' '■U ' * U-' '■fj, 



There are now two new gauge symmetries 



6A^ = V^A, 



A, 



;i.2) 



;i.3) 
;i-4) 



and fixing the gauge A^ = 0, = recovers the original action (1.1). 

We next make the following conformal transformation which serves to unmix the scalar 
and the tensor 



V = Kt'^ jyzr^'^'^'f^di.u- 



The resulting action is 
S = Jd^'x C^=o{h') + 

- 2m^ [h'^.V^A" - h'V^A^) + 2m ^ ^ _ ^ 



\m\h'^,h'^''' - h'^) - Kn" F^,F^^- + '^rn'RA^A^ 

D-1 o R 



m 



D 



(20V^A^ + /i'<; 



where Cm=o{h') is the Lagrangian for a massless graviton and F^^ = d^Ay — dyA^. 
Here we see that for the special value 

(1.7) 

D-2 ' ^ ' 

the dependence on cj) completely cancels out of (1.6). Tracing back through the Stiikelberg 
replacements and field redefinitions, we conclude that the original Lagrangian (1.1) has the 
gauge symmetry 

V^V. + a, (1.8) 

where a is a scalar gauge parameter. This symmetry looks like a tuned combination of 
a linear diffeomorphism with diff parameter ~ d^a, and a linearized Weyl transformation 
with parameter ~ a. This is the partially massless theory [2-9]. Due to the enhanced gauge 
symmetry (1.8), this theory propagates one fewer degree of freedom than the generic massive 
graviton, so in D = 4 it carries four degrees of freedom rather than five. 

To more clearly see what is happening for various values of im? vs. -R, we take the high 
energy limit 

2 

^ 0, i? ^ 0, fixed, (1.9) 

R 

and keep fixed the canonically normalized fields = mA^, = m^0, so as to preserve the 
number of degrees of freedom in the limit. What remains is a decoupled action for a massless 
graviton, massless vector and massless scalar. 



Looking at the signs of the various kinetic terms, we see that the vector field is ghost- 
like unstable^ for < 0. At = 0, where the sign is turning over, only the tensor 
propagates and the action (1.1) has the enhanced linearized diffeomorphism symmetry of 
linearized general relativity 6h^y = V + Vjy^^. 

When m > 0, the vector is stable but the scalar is unstable unless > j^^^R- This 
is the Higuchi bound [3], and when it is saturated the scalar kinetic term disappears and we 
have the enhanced scalar gauge symmetry (1.8) of the partially massless theory. Note that 



^Note that there is no equivalent of the Breitenlohner-Freedman bound [24] for massive gravitons in AdS, 
since we have a ghost-Uke instabihty just as soon as < 0. 
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because of the stability requirement > 0, the partially massless theory is only stable on 
de Sitter space -R > (which happens to be the correct sign of the cosmological constant 
in our universe). In fiat space, ^ = 0, nothing like it exists, and the novel gauge symmetry 
(1.8) merges with the diffemorphisms. See figure 1 for a summary of the various regions in 
the plane. 




Figure 1. Regions of stability for massive gravity on maximally symmetric backgrounds. 

The partially massless theory propagates a particle which lives in an exotic irreducible 
representation of the de Sitter group^ [2-9, 25]. There is no flat space counterpart to this 
representation - in the flat space limit the particle breaks up into a massless graviton and 
a massless photon. They are only seen to be unifled on large distance scales comparable to 
Hubble. 

The partially massless symmetry ties the value of the graviton mass to the curvature 
of the background de Sitter space-time, through the relation (1.7). Said another way, a bare 
cosmological constant is forbidden by the partially massless symmetry. Indeed, at the linear 

^This particle can also be given an electromagnetic- like interpretation [25], complete with duality [26]. 
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level, a bare cosmological constant is \/—g ~ and this is not invariant under the symmetry 
(1.8), bh 7^ 0. In contrast, a large bare cosmological constant is allowed in GR because the 
cosmological constant term h is invariant under linearized diffeomorphisms. At least at the 
linear level then, the partially massless symmetry is a rare example of a symmetry which 
can rule out an arbitrary cosmological constant. 

1.2 The candidate non-linear theory 

Generically, adding higher order interaction terms to the known linear PM Lagrangian will 
destroy the partially massless symmetry and there will be more propagating modes at the 
non-linear level. The question, then, is whether it is possible to have a fully non-linear theory 
of partially massless gravity for a particular choice of interaction terms, which has a scalar 
gauge symmetry like the linear theory. 

The recent developments in massive gravity allow us to revisit this question. The dRGT 
theory of massive gravity [27, 28] provides a non-linear completion of the Fierz-Pauli action 
that is free of the Boulware-Deser ghost [29] at the fully non-linear level [30-39] . The original 
ghost-free dRGT model was proposed with a Minkowski reference metric [27, 28], but this 
was subsequently generalized to a generic reference metric [33, 40]. Just like dRGT massive 
gravity, the PM theory should be free of the Boulware-Deser (BD) ghost instability. Since 
the dRGT models are the unique ghost-free massive gravity theories in D dimensions (at 
least with the conventional Einstein-Hilbert kinetic term, we'll say more about this later), 
then we expect that the PM gravity theory will be a special case of the dRGT massive 
gravity models formulated with a de Sitter reference metric. 

We start by looking for a non-linear PM theory within the class of ghost-free dRGT 
Lagrangians of massive gravity, formulated with respect to a non-dynamical dS reference 
metric g^y, 



Here the dynamical metric is given by g^^y. The symmetric polynomials S'n(M) of a matrix 
M% are 




(1.10) 



(1.11) 
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with 5*0 = 1. The curvature i?, cosmological constant A and Hubble constant H of the 
reference metric (7^1, are related by 

D-2\^ (D- 1){D - 2) 



Not all the /9's are independent: to ensure that g^y = g^y is a solution we enforce 
tadpole cancellation which gives 



D\ V = (D-1)\ V . (1.13) 



(For D = 4, this yields (3o = — (3/3i + 3/32 + /^s)-) We may absorb one further coefficient by 
taking to be the mass squared of the graviton fluctuations around g^^, which is enforced 
by the relation 

Pc — Ji — 2 

(For D = A, this yields jSi + 2(32 + /^a = 8.) Finally, since Sd{,M) = det M the term in (1.10) 
proportional /Sat is a constant and can be ignored, so we set (3^ = 0. In total, the theory has 
D — 2 free parameters, in addition to the mass and cosmological constant. 

We will see in the next section that various lines of evidence (including the earlier 
analysis of [21]) point to a dRGT theory with one particular choice of coefficients which has 
the potential to propagate a partially massless graviton non-linearly. These coefficients are 

/3o = -4(D-l), /32 = -^^, (1.15) 

with all other /3's zero. In addition, the mass should be at the partially massless value 
(1.7). Plugging these into (1.10), and using the relations (1.12) for the background, we have 
our candidate D-dimensional PM Lagrangian: 

CpM = -^V^ [m - S^iv^g)) . (1.16) 
Other than the value of the de Sitter radius, there are no free parameters left in the theory. 



We note that in D = A, we have the property \^—g S2{^/ g ^g) = \f--g 5*2(^/5' ^g) [40]. 
Thus if we promote the PM candidate Lagrangian to a bimetric Lagrangian by introducing 
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an Einstein-Hilbert kinetic term for (7^1,, this new Lagrangian will enjoy a symmetry under 
the interchange of g^.^ ^ g^y. However, we also note that the values of the /3's presented 
here do not coincide with those found in [41, 42], where the authors looked for a partially 
massless theory in a bimetric setup, using different criteria than those used here. 

1.3 PM action in vielbein form 

The non-linear PM action (1.16) can be expressed compactly using vielbeins and differential 
forms'^. We introduce a vielbein 1-form e" and curvature 2-form i?"''[e] for the dynamical 
metric g^J^y^ as well as a vielbein 1-form e"^ and curvature 2-form ^"^[e] = Ae'' for the dS 
reference metric g^y. We invoke the arguments of [38], which allow us to make the following 
replacement, 

y d"" x^g S^^^Y^) -> J - ^ A e'^^^ A ... A e'- . (1.17) 

With this notation, the PM candidate action given by (1.16) can be written in terms of the 
difference of the dynamical curvature and the dS reference curvature, 

SpM = J ea,...aAR^^''' - R'^n Ae'^^A...Ae--. (1.18) 

Let us define the relative spin connection as = Then there are many 

equivalent ways of rewriting the above expression. In particular, we have 

R^b _ j^ab ^^^ab _^a^^^cb (^^^g) 

= m'' + n"^ A n"'' (1.20) 

= 1{V + V)Q^'' . (1.21) 

We use P and V to denote the derivatives which are covariant under local Lorentz transfor- 
mations, for the dynamical metric and reference metric, respectively. 

In this language it is straightforward to demonstrate the existence of the usual linear 
PM symmetry. At lowest order in fields, the above action is invariant under, 

Se" = V {{e'^Y>'d^a) + H'^^a , (1.22) 

with gauge parameter a. 

■^A frame-like formulation of partially massless gravity at linear level was developed in [43]. 
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2 The Evidence 



In this section, we present evidence that the candidate ghost-free dRGT Lagrangian has a 
partially massless symmetry non-linearly. 

2.1 The decoupling limit 

In section 1.1, we saw at the linear level through the Stiikelberg analysis how the longitudinal 
mode disappears from the action at the special value of the mass (1.7). Starting with the 
Lagrangians (1.10), one can perform an analogous analysis at the full non-linear level. The 
Stiikelberg replacement (1.2) must be supplemented by the appropriate non-linear terms [10]. 
Once this is done, the theory is an effective field theory which looks like (1.6) at linear level, 
but now has interaction terms suppressed by various energy scales. The lowest such scale 
(We'll work in D = 4 in this section for ease of displaying the scales) is A3 ~ (Mpim^)^/^ ~ 

(Mpii72)l/3. 

The action can then be greatly simplified by taking the high-energy decoupling limit, 
generalizing (1.9) 

2 

— Tfl 

m^O, R^O, Mp -> 00, ^,A3 fixed. (2.1) 

R 

The decoupling limit focuses in on the leading non-linearities, leaving them intact while 
sending all the others to zero. If the theory propagates a PM graviton non-linearly and has 
only four degrees of freedom, then the mode should disappear from the action completely 
(perhaps after some field redefinitions, generalizing the conformal transformation (1.5) at 
linear level). A necessary condition for this is that the (p mode disappear from the decoupling 
limit Lagrangian. 

This analysis was done in [21] for arbitrary D, and the result is that there is a unique 
choice of the /3„ such that the helicity-zero mode of the massive graviton is absent in the 
decoupling hmit. These coefficients are precisely (1.15), leading to the candidate Lagrangian 
(1.16). 

The analysis of [21] considered only the tensor and scalar modes of the effective theory. 
This is a consistent truncation, but if the partially massless theory exists, the scalars should 
also disappear from the vector-scalar interaction which are present at the scale A3, and 
should in fact disappear completely from the effective theory (possibly after some field re- 
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definitions). In this subsection, we will assume that this is the case, and ask what it implies 
for the resulting effective theory. 

If the scalar has completely disappeared, then the effective field theory only contains 
tensor and vector modes, and the lowest interaction scale that can appear in this theory is 

1 /2 

the one carried by the vector self-interactions, which carries the scale A2 ~ (Mpm) ' ~ 

1 /2 

[MpH] ' . The interaction terms carry one derivative per A, and so they take the form 
J^^A ~ , in terms of the canonically normalized vector field A„ = A2v4„. 

Since there cannot be interactions at energy scales below A2 in a PM theory of gravity, 
one should be able to consider the following limit 

Tfl 

Mpi 00, m,H^O, keeping 77,^2 fixed. (2.2) 

H 

In this limit, one has a flat-space theory propagating a free linearized helicity-2 mode and 
a self-interacting vector degree of freedom, which decouples from the graviton, so this cor- 
responds to a genuine decoupling limit of PM gravity. The decoupled vector Lagrangian is 
symbolically of the form 

^A,dec ~ J^f,,, + + + ■■■ . [2.S) 

At linear level, this Lagrangian propagates only two degrees of freedom since 
respects the Maxwell U{1) gauge symmetry. Non-linearly, the vector field decoupling limit 
Lagrangian should continue to propagate only two DOF. This implies that the interactions 
should remain gauge invariant under some gauge symmetry which starts as the Maxwell 
symmetry and then possibly gets corrections at higher order in A, 

6A^ = d^A + --- . (2.4) 

To see whether such a symmetry can ever exist non-linearly, let us start with a massive 
graviton g^^, on an (A)dS background g^j,^. To make the argument as general as possible, we 
work in what follows perturbatively with an arbitrary mass term (although we will see soon 
how the mass term is constrained to the coefficients of (1.16), at least up to cubic order in 
perturbations). 
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Thus our starting point is a Lagrangian containing the Einstein-Hilbert kinetic term 
with cosmological constant A, and an arbitrary mass term 



■£ — -^^EH + 1^ 



2 



-g{R\g\ - 2A) - —^g {\h^\ - \hf + c^\h^\ + c^\h\h\ + C3[/i] 



2 

+ d2\h''\\h\ + 4[/iT + d^\h\h\^ + ci5[/i]' + ■■■)]' (2-5) 

where the arbitrary mass term is expressed in powers of h^j_y = g^y — g^y, and indices on h 
are raised and lowered with g^^. The mass term is the Fiertz-Pauh term at quadratic order, 
followed at higher order by every possible contraction of h^i, with an arbitrary coefficient in 
from of each. 

To derive the contributions to the Lagrangian (2.3) in the decoupling limit we make 
the replacement [10] 

Ku ^ d^Ay + dyA^ - d^APdyAp, (2.6) 

in the mass term, and take the metric and all covariant derivatives to be flat. We can do 
this because by hypothesis all the couplings to h disappear in the decoupling limit (2.2), as 
do any further corrections due to the (A)dS curvature of the reference metric. 

We expand the Lagrangian (2.3) in powers of A, 

C-A = ^{2) + ^{3) + ^{4) + ■ ■ ■ , (2.7) 



'PI' 

and look for a scalar gauge transformation with parameter a, in powers of A, 

SA, = L,a, L, = Lf + L« + + ■ • • . (2.8) 
If the Lagrangian is gauge invariant, we should have the Bianchi identity 

4^^0. ,2,9) 



At lowest order, quadratic in A, we have £(2) = — |-F^^, which is invariant under 



At cubic order in A, the Lagrangian (2.5) with the replacement (2.6) gives 
1 



(3) 



-2 + 3ci) dyAp dpA" d^A^ + ci dyA^ O^Ap dpA" + 4c3 {d ■ Af (2. 10) 
+ 2c2 {d ■ A) d'A^ d^Ay + 2(1 + C2) {d ■ A) d^A" d^A, . 
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Next, we write the most general transformation L^f^\ first order in A, with up to three 
derivatives (though only up to two actually appear, because we can't contract three), with 
an arbitrary coefficient in front of each term, 

L«a = Cl A^a + C2 A^^Ua + C3 A'^df.d^a + C4 d^A^d^a 

+C5 d'A^ d,a + Ce (9 ■ A) d^a + Cj d.AA" a + Cs DA^a . (2.11) 



At third order in A, the identity (2.9) reads 

f (0) ^^(3) , f (1) ^^(2) 



6A^ ^ SA^ 
Demanding this hold fixes the cubic mass term to the values, 

Cl = 1, C2 = Cs = ^ , (2.13) 

which correctly matches the cubic mass terms of our candidate theory (1.16). The gauge 
transformation is fixed to be 

L«« = ^F.^d'^a + d, [C^A^d'a + C^{d ■ A)a] . (2.14) 

The unfixed coefficient Cq corresponds to the freedom to redefine the gauge parameter a — 
[d ■ A) a and C3 corresponds to the freedom to redefine the gauge parameter a — )-~ Aj^d^a. 

While we displayed the results of this subsection for D = 4, this analysis does not 
depend on the number of dimensions and the coefficients (2.13) are the same for any D, 
and match those of (1.16) for any D. 

In D = 4, the effective ffeld theory of PM gravity at distances shorter than Hubble is 
a theory of a massless tensor and a vector. If we take the graviton mass to be the observed 
value of Hubble today, and the Planck mass to be what is observed, then the cutoff scale 
comes out to 

A2 ~ neutrino mass scale. (2.15) 

This cutoff is parametrically higher than the cutoff A3 generically present in dRGT massive 
gravity, because the worst offending scalar interactions are absent. It's worth noting that 
this is the sort of effective theory one would obtain by integrating out all massive particles 
from the Standard Model; the only massless degrees of freedom are the photon and graviton, 
and the cutoff is at the neutrino scale, where the first massive particle comes in. 
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2.2 Non-linear gauge symmetry in mini-superspace 



The candidate action (1.16) propagates a partially massless mode at linear level, the scalar 
mode is absent in the scalar tensor sector of the A3 decoupling limit, and it has the required 
gauge symmetry at cubic order in the A2 decoupling limit. The key question we would like to 
answer is whether or not the PM symmetry exists in the full theory, with all non-linearities 
and beyond the decoupling limit. While this is a difficult question to address at the fully 
non-linear level for an arbitrary metric, we can simplify the analysis by adopting a particular 
ansatz for the dynamical metric and checking for the existence of the PM symmetry there. 

Here we study the full non-linear theory in the simplified case of mini-superspace, i.e. 
for an FRW ansatz. Consider the PM Lagrangian in the form (1.16), only with an arbitrary 
coefficient A in front of the mass term, so that we can see what happens as A — iJ^, 

S = J d^x ^ (^R[g] - 2XS,{^/^g)j . (2.16) 

We plug in a flat FRW ansatz in the flat slicing for the reference dS metric, with a dynamical 
scale factor a{t) and lapse N{t) (note that there is no diffeomorphism symmetry to set the 
lapse to unity), 

g^^ = -N{tydt'^ + a{tydx\ g^^ = -dt^ + e^"'dx^ . (2.17) 
Up to a total derivative, the action becomes 



S = M^f^ I dt 



N 



2A— — a^-'e^* a + ——Ne 



(2.18) 

The lapse N{t) appears algebraically, so we may eliminate it through its equation of motion, 

Ar=4^e-^*d. (2.19) 



Plugging back in the Lagrangian, we flnd up to a total derivative, 

S = J d^x {D - l)yAe^*a^-2(i/ - VX) . (2.20) 

For the generic case A 7^ H^, the a equation of motion tells us a = and the Lagrangian 
is inconsistent (which is an instance of dRGT massive gravity having no flat FRW solutions 
[44]). But for the special case 

A = H\ (2.21) 
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the Lagrangian becomes empty, meaning there is a gauge symmetry with gauge parameter 
e{t) that allows us to set the scale factor to anything we please, 

5a = e, 6N=^e-^'e, (2.22) 

In this theory, cosmology is pure gauge: any arbitrary function is a solution for the scale 
factor, with the lapse determined by (2.19). The transformation (2.22) should be the mini- 
superspace form of the full non-linear PM gauge symmetry. 

In terms of the metric, this implies that the PM Lagrangian is invariant under the 
following transformation, 

= diag(-2iv4e-^*e, 2ae, ■ ■ ■ ) • (2.23) 
II 

The transformation here is only first order in derivatives. If we perform the change of 
variables e = ^e^'^{—Ha + H^a) we find that this transformation is entirely consistent with 
the linear PM symmetry (1.8) with gauge parameter a{t). 

2.3 Perturbations in FRW 

Following [45], let us consider cosmological perturbations around FRW for dRGT massive 
gravity on dS. Requiring the coefficient of the kinetic term of the helicity-zero mode to be 
positive is what gives rise to the generalization of the Higuchi bound [3] . Specializing to four 
dimensions D = 4, this was derived in [45] to be (we convert the a„ coefficients used there 
to the /3„ coefficients used here) 

m\H) = ~[^i + 2/3,^ + /33^) > 2H^ , (2.24) 

where now H is the Hubble parameter of the dynamical metric, and Hq that for the reference 
metric. The dynamical quantity im?{H) is the actual mass of the graviton for the fluctuations 
when the dynamic metric and the background metric are not the same. 

If the PM symmetry exists then the kinetic term for the helicity-zero mode should 
vanish altogether. This corresponds to saturating the bound. However the bound must be 
saturated for all values of H since this bound is applicable to any spatially flat FRW geometry 
[45] . The only way this is possible is if the coefficient of each order of the polynomial in H 
balances on each side of the equation. That is in order to satisfy 

m\H) = ~U + 2P2§- + /^sfj) = 2H^ (2.25) 



- 16 - 



for all H we must satisfy 

/3i = (2.26) 

(32 = 8H^/m^ (2.27) 

/33 = (2.28) 

Using the fact that the PM condition requires = 2Hq, we find (32 = 4. Along with 
the tadpole cancellation condition (1.13), this gives (3o = —12. These are precisely the 
same coefficients as those of our candidate theory (1.15) in D = 4. I.e., they are the same 
coefficients as the those found in the decoupling limit analysis of [21] (expressed here in 
terms of (3n)- 

In other words, adopting the usual Higuchi bound condition = 2Hq, along with 
the same coefficients (1.15) of the dRGT model found in the decoupling limit analysis, 
guarantees that the kinetic term for the helicity-zero mode vanishes automatically for any 
FRW geometry. That the cancellation takes place for this specific choice of parameters 
demonstrates a nontrivial consistency check. This is an independent confirmation of the PM 
candidate Lagrangian, performed away from the decoupling limit. 

2.4 The full theory at cubic order 

In this section we will perform a brute force systematic analysis, away from the decoupling 
limit and for an arbitrary dynamical metric, so as to determine order by order in non-linearity 
whether or not a PM theory of gravity exists. Our starting point is a Lagrangian containing 
the Einstein-Hilbert kinetic term with cosmological constant A, and an arbitrary mass term 
written in powers of h^i, = g^y — g^y. Indices on h^y are raised and lowered with the reference 
metric g^^y, 

m 



'PI 



-g{R{g]-2K)-'^^g{C^l^ + Ct^ + Ct^ + ...) , (2.29) 



with 



C^^) = h,[h^] + h2[h]\ (2.30) 
= ci[/i3] + C2[h\h] + C3[/^]^ (2.31) 
C^^) = di[h^] + d2[h^][h] + + d,[h''][h]^ + d,[h]^ , (2.32) 



-17- 



Here we do not restrict ourselves to the dRGT interactions - the mass term is simply every 
possible contraction of hfj,^, with an arbitrary coefficient in from of each. Again, the curvature 
of the reference metric R is related to the cosmological constant A by (1.12). 

If a PM theory of gravity exists, it must have a scalar gauge symmetry, and be invariant 
under a transformation of the form 



Sh 



L ni/Oi • 



(2.33) 



where L^,j is some operator and a is the gauge parameter. Gauge invariance then gives us 
the Bianchi identity 



5C 



6h 



111/ 



1 



liu\ 



6h 



0. 



(2.34) 



The operator L can be expanded in powers of h, 

L = + + ■ ■ • . 



(2.35) 



In what follows we will attempt to determine the coefficients of the mass term and of the 
terms appearing in the perturbative expansion of L, in order to solve (2.34) order-by-order''. 



2.4.1 Quadratic order 

At lowest order, the statement for gauge invariance (2.34) reads 



-g (G'^'^ + A<7ni(i) + 



6h 







(2.36) 



where here and in what follows, the subscript designates the expansion to n^^ order in 
hfiu- We restrict ourselves to gauge symmetries with at most two derivatives. The most 
general transformation which does not involve any power of the metric perturbation then 
has three possible terms: 



(2.37) 



where V represents the covariant derivative with respect to the background dS metric, and 
we have put an arbitrary coefficient in front of each possible term. 

^The is the unitary gauge version of a similar analysis done in Ref. [22], and we find equivalent conclusions. 
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Plugging the transformation (2.37) into (2.36), we find that the coefficients of the 
mass terms and of the transformation law get fixed. To see this explicitly we can first fix 
coefficients starting from the highest derivatives: requiring terms of the form VVVV/i in 
(2.36) to vanish fixes = 0. Then we work our way down, requiring terms of the form 
VV/i to vanish, and then zero derivative terms h. This procedure fixes 

2A „ 2A 



h = -62 



-Bi, B-i 



0. 



(2.38) 



D-r ^ {D-i){D-2y 

(Notice that there is also a solution with 61 = 62 = -B2 = -B3 = 0, which just corresponds to 
massless GR with the scalar part of diffeomorphism invar iance.) Scaling i?i = 1 by absorbing 
it into the gauge parameter, we find the 'semi-conformal' transformation of (2.37): 

2A 



^^-=^'^^^"+(D-l)(D-2) 



(2.39) 



We find that the symmetry fixes the graviton mass to the partially massless value 

2A 



m 



D 



(2.40) 



Finally, the symmetry fixes the Fierz-Pauli structure: bi = —62 = 1- The is a sign that the 
partially massless symmetry knows about the absence of ghosts. 

2.4.2 Cubic order 



At cubic order, the statement for gauge invariance (2.34) reads 



^(0) 

fjLV 



(3)' 



= 0. 

(2.41) 



For I/^^-*, the most general form up to two derivatives has 18 terms (again all indices 
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moved with which we write with arbitrary coefficients: 

= 

+C4 V^^K)^Vxa + V(^/iV,)a + Cg Vxh^u^^a + Cr Vx\^^V,)a 
+C8 V^/i a + Va V(^/i^)^ a + Cio □ V « 
+C11 g^^^hUa + C12 ^^^/i^'VAV^a 
+C13 ^^.V^/i'^Vao + Ci4 ^^.VA/iV^a 
+C15 g^^u^h a + Va V<,/i^" a 

+C17 + C18 g^,uha. (2.42) 

These are organized as follows: We start with terms containing two derivatives. First are the 
terms which are not proportional to g^p. The first line has the terms with both derivatives 
on a, the second line has the terms with one derivative on a, and the third line has the terms 
with no derivatives on a. Next we write the terms proportional to g^^. The fourth line has 
the terms with both derivatives on a, the fifth line has the terms with one derivative on a, 
and the sixth line has the terms with no derivatives on a. 

We substitute this general expression into (2.41) and fix the coefficients Ci for the gauge 
transformation as well as the coefficients q in the mass term by demanding that the equality 
(2.41) holds. Requiring the various terms of the form /iVVVV/i (6 possible contractions) 
to vanish fixes Ci = Cn = C12 = 0. Then we work our way down, requiring terms of 
the form VhVVVh to vanish (11 possible contractions), then terms of the form VVhVVh 
(13 possible contractions), then the two derivative terms hVVh (5 possible contractions) 
and VhVh (5 possible contractions), and finally the zero derivative terms h"^ (2 different 
contractions), fixing coefficients along the way. We have 21 free coefficients to kill these 42 
terms, so we need 21 miracles. 

The cubic mass term coefficients are fixed to the values 
1 



2 



^g{R[g] - 2A) - ^^=1 ^[h'] - [h]' - [h'] + ^[/.^][/.] - \[hf + 



(2.43) 



with fixed to (2.40). The coefficients of the mass term determined by this procedure 
agree with those found in an expansion of the PM candidate Lagrangian (1.16). 
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The gauge transformation is fixed to be 

1 - ^ A D-Q 



+C3 



- - 2A _ 

V^V^ {ha) + p_ !)(/}_ 2) ^^'' 



(2.44) 



There is one unfixed parameter C3, which is just the hnear transformation again with pa- 
rameter ha, so by re-defining the gauge parameter we can set C3 = 0. 

All the needed miracles occur only D = 4, in other dimensions we are two short^. The 
zero derivative terms in (2.41) are not canceled, except in D = 4, 

-,„,i£<^ 



Note that the obstruction for D 7^ 4 vanishes in the decoupling limit, and curiously, 
it has the Fierz-Pauli structure. This seems to imply that the partially massless theory as 
we have it can only exist in D = 4 dimensions. As has been noted, D = 4 is special in that 
partially massless theories are conformally invariant there [8, 22]. Also, this is the dimension 
in which the partially massless candidate has a Z2 symmetry between the dynamical and 
background metrics. 

2.4.3 The cubic gauge transformation 

At quadratic order, the derivative part of the transformation had the form of a diffeomor- 
phism, i.e. we could write 5/i|fj|deriv = V^^i-^^ + Vjy^^°\ where ^^"^ = ^V^a. We can ask 
if the derivative part of the cubic transformation has the same property, that is, if we can 
write the non-linear transformation determined thus far as an expansion of 

Sh^uUeriv = V^^i. + Vu^^ , (2.46) 

for some vector 

To expand, we write the covariant derivative V in terms of the background covariant 
derivative V and an associated connection C^^^ 

V^e. = V^e.-C^V^A, (2.47) 

^Curiously, the prediction D — 4 seems to agree with observations. 
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= b^'^iSf^gu. + V,(7^« - V^g^,) . (2.48) 

Expanding (2.46) in h^y gives 

Sh^il Ideriv = V,e^ + V.e?) - 2CfJif , (2.49) 

for the transformation beyond leading order. 

As determined by the analysis above, we have 

Idcriv = |/iA(MV.) V^a + I Va VV^a - | V(^/i,)AV^a , (2.50) 

for the derivative part of the cubic transformation. Equating (2.50) with (2.49), we can 
solve: 

^^ = iVV'«- (2-51) 

Thus, indeed, to cubic order the derivative part of the transformation has the form of a 
diffeomorphism with parameter 

^z. = iV^a + i VV^a + Oih^) . (2.52) 
3 Obstructions to Non-Linear PM Gravity 

In this section we discuss obstructions to a general non-linear PM symmetry from several 
different angles. In particular, we show that the candidate PM theory (1.16) does not have 
the required gauge symmetry to all orders. 

3.1 Anisotropic obstruction 

We now consider an extension of the mini-superspace analysis of Section 2.2, where we still 
restrict ourselves to an homogeneous ansatz but allow for anisotropics in the physical metric, 
so that the metric (2.17) gets replaced by 

D~l 

g,u = -N{tfdt' + ^ a'itfdx'' , (3.1) 

i=l 

while the reference metric Qf^i, remains dS. We will show below that, while the isotropic ansatz 
had an evident gauge symmetry, this is not the case as soon as one allows for anisotropics. 
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This represents a first hint towards the absence of a PM symmetry in the fully general and 
non-linear case. 

Let us focus on the case that _D = 4 so that our candidate action is given by 

S=^ jd^x^—g (R[g] - 2XS,{y^g)) . (3.2) 

Using our anisotropic ansatz (3.1), and after integrations by parts to remove higher deriva- 
tives, the PM action in four dimensions is 



AT ''J^ V ^ 

We have defined the fully symmetric constant arrays Sl^\ Slj\ S^jl so that 



(3.3) 



Sl^'a' = a^ + a^ + a^ (3.4) 



Slfa'a^ = a^a^ + a^a^ + a^a^ (3.5) 
Sljla'd^d'' = a^d^a^ + a^d^d^ + a^d^d^ . (3.6) 

We also define the shorthand S^^^ = ^a', ^(2) ^ S^fa'a^, and S^^^ = Sl^la'd^dK 

Similarly as in the isotropic case, the lapse is an auxiliary variable and we may solve 

for it, 



-Ht 



^-vlVsoy (3.7) 

Plugging this constraint back into the Lagrangian, we get 

C = -M|ie^*yA [2VS'(i)5'(3) + ^S^^A . (3.8) 

We want to know whether this Lagrangian has a gauge symmetry, so we proceed with 
a Hamiltonian analysis. The canonical momenta are 



vr. = g = -2M|,e^*yA^ fi^S«^«^ • (3-9) 

Defining the 'phase space metric' Gij as Gij = S'^'^^ S^^^la^ , and G^^ its inverse, we see that 
the canonical variables satisfy TCiG^^TTj = AMp^e'^^^X , so we have the primary constraint^ 

C = ^TTiG^^TTj - 2Ml^e"'V\ . (3. 10) 



^We have square rooted both sides because this is the constraint that is regular in the isotropic case, and 
we have chosen the sign because the canonical momentum is negative (3.9). 
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The Hamiltonian is given by 



n = 'K,a} - C = Ml^e^'\S^^\ (3.11) 

and does not depend on the canonical momenta. 

If the primary constraint (3.10) is associated with a sole gauge symmetry, it should be 
the only constraint (in which case it is first class and generates the gauge symmetry). Thus 
there should be no secondary constraint. To look for a secondary constraint, following Dirac, 
we check the primary constraint's conservation in time, 

f^^^ ^C,n} = f - f = -2Mg.e«-VA (i/ + ^'^'f"'] . (3.12) 

dt dt ^ ' ^ dt dni da' \^ y/TCiG^^nm J 

In the generic anisotroic case, the right hand side does not vanish on the constraint surface, 
so we have a secondary constraint and hence not the single gauge invariance we are seeking. 



In the isotropic case, G^^^jSl^ a!^ ~ a ~ ~\/'^eG^"^Tim (remembering tt < 0), and the 
right hand side vanishes when positing A = H^. However we see now that this this gauge 
symmetry is lost as soon as we are dealing with a more general metric, even when A = H^. 
This is the first hint that the symmetry found in the previous section, either fully non-linearly 
in the minisuperspace, or at linear level does not propagate in full generality in the most 
naive version of the theory. 

3.2 Quartic obstruction from vectors in the decoupling limit 

We now move onto a second way to see that PM gravity fails to be realized at the non-linear 
level. For this, we push the analysis of the vector modes in Section 2.1 to quartic order in 
the fields. We will see that the required gauge invariance of the vector modes cannot be 
realized at quartic order. 

At quartic order in the vector field, the requirement (2.9) of gauge invariance for the 
decoupled vector mode Lagrangian reads 

^r^ + ^l^^^ + ^r^ = 0- (3-13) 
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The most general transformation Lfj, which is at most second order in derivatives and 
quadratic in the vector field is, 

L(^)a = Di A^d^a + D2 A^d^a + d^A^ A'' a 

+D4 d,A^ A" a + Z^5 (^^) «• (3-14) 

Plugging this into (3.13), and using and using the lower order relations from Section 
2.1, we find that no matter what the coefficients Di, ■ ■ ■ , of the transformation are, and 
no matter what choice of potential [i.e. choice of coefficients (ii, ■ ■ ■ ,^5) allows (3.13) to 
be satisfied^ This remains true if we allow for up to 3rd derivatives in the ansatz (3.14) (as 
might be expected from counting derivatives in the decoupling limit. 

We therefore conclude there is an obstruction to preserving only 4 DOF at the quartic 
level when considering a generic mass term in arbitrary dimensions. It is worth pointing out 
that this obstruction relies on the assumptions that a) the Einstein-Hilbert term receives no 
corrections, and b) that the gauge transformation is at most 3rd order in derivatives when 
acting on A. It is possible that the obstruction could be relaxed by allowing more derivatives, 
however in that case we would also have more conditions to satisfy and it could be that the 
transformation would then involve an infinite number of derivatives and would thus likely 
be non-local^. 

''In particular, if we look at the two tensor structures A" d ^A^^d pAy and d'' A*^^^lAP^^pAy coming 
from (3.13), we find their coefficients depend only on d\ and di in the following way 

(I _ 2di - dt^A'^d^^APUdpAy + Q - 2di - ^-d^ d'' A^dpAPUdpAy, (3.15) 

and so both cannot be made to vanish simultaneously. 

®Note that this analysis of this section also serves to rule out the possibility that the flat-reference metric 

"minimal model" (in the language of [46]), for which the decoupling limit interactions vanish, has the higher 

interaction scale A2 (in D — K). This is because if the theory had this higher cutoff, the vector self- interactions 

in the A2 decoupling limit would take the same form as they do in our analysis here, and since they cannot 

have a gauge symmetry non-linearly, they propagate more than two degrees of freedom non-linearly. Then 

the decoupling limit theory (remembering the scalar kinetic term is still present) has more than the 5 degrees 

of freedom required of the ghost-free massive graviton. This is consistent with the amplitude analysis of [47] 

which shows A3 as the cutoff, which itself implies that the vector-scalar interactions in the decoupling limit 

of the minimal model cannot vanish. 
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3.3 The full theory at quartic order in D = 4 



Despite the encouraging evidence of section 2 for the existence of a PM theory of gravity 
propagating only four degrees of freedom in four dimensions, the previous subsections aheady 
provided a few arguments obstructing the existence of a fully fledged PM theory of gravity. 

Here, we report on our attempt to push the brute force calculation of Section 2.4 to 
fourth order in the flelds in D = 4 (we were obstructed already at cubic order when D A). 
We flnd that there is an obstruction, and the theory fails to be gauge invariant at this order. 

At quartic order, the statement for gauge invariance (2.34) reads 



1 

'2^ 



1 

'2^ 



-9 {G^" 



5C 



(4) 



-9 (G 



5C 



(3) 



fii/ 



6C 



(2) 



'(1) + 5h 



0. 



(3.16) 



For L^'^\ the most general form up to two derivatives has 72 terms (which we will not attempt 
to display here), with the same number of free coefficients. In addition, we keep the quartic 
mass term general, which brings 5 free coefficients di, - ■ ■ ,(1^. We plug this into (3.16), using 
the quadratic and cubic order results (2.39) and (2.44) for the lowest order transformations, 
and we demand that (3.16) holds. We fix coefficients starting from the highest derivatives, 
fixing coefficients requiring the various terms of the form /i^VVVV/i to vanish, then the 
terms hVhV'V'Vh, then V/iV/iVV/i. All these four derivative terms can be eliminated. 
Next we start on the two derivative terms, starting with those where both derivatives act on 
the same h, i.e. those of the form /i^VV/i, which can also be eliminated. Then we come to 
the two derivative terms of the form hVhS/h, and it is here where an obstruction appears: 
these terms cannot be canceled without setting A = 0. 

In the process, however, the coefficients cii, ■ ■ ■ , (is in the quartic mass term in (2.5) 
are all uniquely determined. 



C 



m 

1 



^?(i?M-2A) 



7 



11, 



m 
T 

212 



1 



32 
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with fixed to (2.40). These mass term coefficients match precisely those coming from 
expanding the candidate theory (1.16) to quartic order (even though we started with the 
most generic non-derivative interactions). Nevertheless, it seems impossible to preserve a 
symmetry at quartic level and hence to ensure that the helicity-0 mode of massive gravity 
gets consistently removed, inhibiting the existence of a partially massless theory of gravity. 
The only assumptions we made are that the kinetic term is that of Einstein-Hilbert, and 
the gauge symmetry has at most two derivatives. This conclusion is mentioned without 
calculation in [22], and reached as well with different methods in [20]. In the case of D = 3, 
new massive gravity [48] on a de Sitter background also becomes a situation where the PM 
symmetry is present at linear level but does not extend to all orders [49]. 

4 Outlook 

Unfortunately, despite the suggestive evidence presented in the first half of this paper, a 
brute force calculation demonstrates that, as formulated, the PM theory doesn't exist at the 
non-linear level. Let us emphasize that while much of our analysis focused on a particular 
candidate theory (i.e., a specific choice of dRGT massive gravity), our "no-go" result is more 
general. In our perturbative analysis we placed no assumptions on the form of the non- 
derivative interactions and we found that no gauge symmetry can exist beyond cubic order, 
in any number of space-time dimensions. 

There are still several ways in which the no-go result may be avoided. In dRGT the 
kinetic term is that of Einstein-Hilbert, and even in our more general perturbative analysis 
we assume an Einstein-Hilbert kinetic term as well. It could be that there are non-canonical 
or higher derivative kinetic terms, analogous to the Lovelock terms, which are still ghost- 
free and allow for partially massless interactions. Another possibility is that additional 
degrees of freedom are needed in order to have consistent interaction. It could be that one 
of the bimetric or multi-metric versions of dRGT theory [38, 40] has a gauge symmetry and 
contains a partially massless graviton. Candidate bimetric theories are proposed in [41, 42]. 
In this case, the PM bimetric theory should describe a massless spin-2 field interacting with 
a partially massless spin-2 field (or additional massive fields, in the multi-metric case). 

A theory which comes close to realizing a non-linear partially massless graviton is 
conformal gravity in D = A, whose action is the square of the Weyl tensor, as studied 
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recently in [50]. Around a de Sitter solution, the theory propagates a massless graviton and 
a partially massless graviton [51], but unfortunately precisely one is always a ghost. Any 
attempt to truncate out this unwanted degree of freedom fails at quartic order in the fields 
[50, 52]. Nevertheless, the full theory has six degrees of freedom and is Weyl invariant, which 
acts as the partially massless symmetry (1.8) at linear level. Thus this theory serves as one 
example where it is possible to extend the PM gauge symmetry to all orders. 

Finally, if PM is to be a phenomenologically viable theory of gravity, the issue of 
how the graviton couples to matter must be addressed. Requiring that the matter sector 
maintains the partially massless symmetry will certainly put strong constraints on the form 
of the coupling. For instance, at linear level the symmetry (1.8) requires that the stress 
tensor satisfy 



So the stress tensor does not have to be conserved, but its non-conservation is suppressed 
by Hubble scales. Nonetheless, the possibility of addressing the old cosmological constant 
problem through the partially massless symmetry remains a compelling reason to continue 
the pursuit of PM gravity and its possible couplings. While ultimately unsuccessful, the 
candidate theory of PM gravity studied here has enough promising features to serve as a 
starting point for future investigations. 
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